Non-Gaussianity of quantum fields during inflation 



- 1—1 

X 



Kazuya KoyamE0 
Institute of Cosmology & Gravitation, Dennis Sciama Building, 
University of Portsmouth, Portsmouth P01 3FX, UK 

In this review, we discuss how non-Gaussianity of cosmological perturbations arises from inflation. 
After introducing the in-in formalism to calculate the n-point correlation function of quantum fields, 
we present the computation of the bispectrum of the curvature perturbation generated in general 
single field inflation models. The shapes of the bispectrum are compared with the local-type non- 
Gaussianity that arises from non-linear dynamics on super-horizon scales. 
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5-H 1 I. INTRODUCTION 

' There is currently a great deal of interest in the statistical properties of primordial perturbations from 

inflation, because measurements of any non-Gaussianity will improve by about an order of magnitude over 
the next few years, for example with the Planck This will provide a key way to discriminate between the 
many models of inflation. Although single field models of slow-roll inflation typically generate a small level 
of non-Gaussianity 0, Ej[ , there may be an observable level generated in many alternative models of early 
universe [3-l58j. 

O |' We are interested in the primordial curvature perturbation on uniform density hypersurfaces, £, on large 

scales, which is directly related to temperature anisotropics in Cosmic Microwave background (CMB). The 
, ■ power spectrum of C is defined as 

(N (C(ki)C(k 2 )) = (2w) 3 6^(k 1 + k 2 )P c (fcx), (1) 

> 

where C(k) is a Fourier component of £(i, x l ). If ( obeys Gaussian statistics, the power spectrum determines 
all statistical quantities. However, if the distribution function of £ deviates from Gaussian statistics, we need 
to specify higher order statistics. The first non-trivial statistics is the bispectrum defined by 

CS : (C(ki)C(k 2 )C(k 2 )> = (27r)V 3 )(k 1+ k 2 +k3)B c (fci,fc 2 ,fc3). (2) 

o . 

Currently observations of the CMB have concentrated on constraining the 3-point function (bispectrum) 

The most instructive way to understand how non-vanishing bispectrum appears from non-linearities is to 
use the so-called delta-N formalism (62l - [65| . See (6^ for a concise review of the delta-N formalism. This 
is based on the separate universe approach [67l469l |. This considers each super-Hubble scale patch to be 
evolving like a separate Friedman-Robertson- Walker universe which is locally homogeneous. By patching 
these regions together we can track the evolution of the curvature perturbation on large scales just by using 
background quantities. 

The number of e-foldings, N, given by 



N = J H(t)dt, (3) 

is evaluated from an initial flat hypersurface to a final uniform-density hypersurface. The perturbation 
in the number of e-foldings, 5N, is the difference between the curvature perturbations on the initial and 
final hypersurfaces. We wish to calculate primordial perturbations, hence we pick a final uniform density 
hypersurface to be at a fixed time during the standard radiation dominated era, for example during primordial 
nucleosynthesis. The initial time is arbitrary provided it is after the Hubble exit time of all relevant scales. 
It is often convenient to pick this time to be shortly after Hubble-exit. 
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Let us consider a model where a scalar field <f> determines the expansion of the Universe. The scalar 
field acquires quantum fluctuations 8<fi under horizon scales which cause fluctuations in e-foldings in each 
super-Hubble scale patch. Then £ can be written as 

C{t,x l ) = N^t,U)Scj)(U,x l ) + ijV,^(t,t») (Scb{U,x 1 ) 2 - (5cj)(U,x 1 ) 2 )) + (4) 

where t* denotes the horizon crossing time and = dN/d(j> and Nm — d 2 N/d<j) 2 . Using Eq. (QJ, the 
bispectrum of the curvature perturbation is calculated as 

(C(k 1 )C(k 2 )C(k 3 )) = Arf^kO^)^ (5) 
where [6<j) * 5<p] denotes a convolution 

[<ty*<ty](k) = J -03«50(q)^(k-q). (6) 

A diagrammatic approach to compute the higher order correlation function is developed in Ref. [70| . 

There are in general two contributions to the bispectrum of the curvature perturbations. One 
arises from the bispectrum of quantum fluctuations of the scalar field generated under horizon scales, 
(<50(ki)<5</>(k2)<50(k3)). The other is coming from the non-linear evolution of the scalar field on super-horizon 
scales determined by the second derivative of N with respect to the field, N^$. The latter contribution exists 
even if the field perturbations at the horizon crossing are Gaussian. This contribution is often called the 
local type non-Gaussianity as this arises from a non-linear local relation between the curvature perturbation 
£ and the field perturbations 5<j) in real space. On the other hand, the contribution from the non-linearity 
of quantum fields depends on non-linear interactions under horizon scales and its fc-dependence is in general 
very different from the local type non-Gaussianity. In this review, we derive this contribution in several 
inflation models and compare its shape with the local type non-Gaussianity. 

The structure of this review is as follows. In section II, we review the in-in formalism to calculate the n- 
point function of quantum fields. In section III, the effective action for higher order perturbations is derived 
in two gauges that are necessary to compute the interaction Hamiltonian at third order. We calculate the 
bispectrum of the curvature perturbation in slow-roll inflation and k-inflation. In section V, we compare 
the shape of the bispectrum in k-inflation models with that in the local-type non-Gaussianity. Section VI is 
devoted to conclusions. 



II. QUANTUM CORRELATIONS IN THE IN-IN FORMALISM 

In this section, we review the in-in formalism to calculate the n-point function of quantum fields by 
following Ref. 71]. Consider a general Hamiltonian system, with canonical variables </>(x, t) and conjugates 
7r(x, t) satisfying the commutation relations 



(x,t),7r(y,t)J = i<5 3 (x-y) , 
and the equations of motion 

0(x,t)=i[ff[0(t),7T(t)],0(x,t) 



<Kx,t),</>(y,t) = 7r ( x ,*) J 7r (y,*) 



= o, 



7r(x,i) =i ff[0(t),7r(t)],7r(x,i) 



(7) 



(8) 



The Hamiltonian H is a functional of the </>(x, t) and 7r(x, t) at fixed time t, which according to Eq. ((SJ is 
independent of the time at which these variables are evaluated. 

We assume the existence of a time-dependent classical number solution <fo(x, t), 7To(x, t), satisfying the 
classical equations of motion and we expand around this solution, writing 



(x,t) = </>o(i) +<ty(x,i) , 



7r(x, t) = 1To(t) + (*>7r(x, t) . 



(9) 



Since classical numbers commute with everything, the fluctuations satisfy the same commutation rules ([7]) 
as the total variables. 
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Now, although H generates the time-dependence of 0(x, t) and 7r(x, t), it is H rather than H that generates 
the time dependence of <5</>(x, t) and 8tt(x, t) where H[S(f>(t),S7r(t);t] is the sum of all terms in H[(f>o(t) + 
5<fi(t),iro(t) + 5ir(t)} of second and higher order in the <5</>(x, t) and/or <5-7r(x, t): 



c,t)=i H[S(f>(x,t),8Tr(x,t);t],S(f)(x,t) 



^7r(x, t) = i H [50(x, t), <5-7r(x, t); t], <5-7r(x, i) 



(10) 



This then is our prescription for constructing the time-dependent Hamiltonian H that governs the time- 
dependence of the fluctuations: expand the original Hamiltonian H in powers of fluctuations 8(f) and Sir, and 
throw away the terms of zeroth and first order in these fluctuations. It is this construction that gives H an 
explicit dependence on time. It follows from Eq. (|10[) that the fluctuations at time t can be expressed in 
terms of the same operators at some very early time to through a unitary transformation 



8(f>(t) = U-\tM)S4>%)U{t,t Q ) , 6%(t) = ir 1 (t,i )a7r(to)l7(t,io) 
where U(t, to) is defined by the differential equation 



(11) 



-U(t, t ) = -i H[5<f>(t),Sw(t); t] U(t, to), 



(12) 



and the initial condition U(to,to) = 1. In the application that concerns us in cosmology, we can take 
to = — oo, by which we mean any time early enough so that the wavelengths of interest are deep inside the 
horizon. From now on, we will omit the x dependence of 5(j) and 5ir to simplify the notation. 

To calculate U(t,to), we now further decompose H into a kinematic term Hq that is quadratic in the 
fluctuations, and an interaction term Hi: 



H[6<j>(t),6w(t);t] = H [6<f>(t),6w(t);t] + H^S^t), 5ir(t); t] 



(13) 



and we seek to calculate U as a power series in Hi. To this end, we introduce an "interaction picture": we 
define fluctuation operators 8(f) 1 (t) and Sir 1 (t) whose time dependence is generated by the quadratic part of 
the Hamiltonian: 



8^ (t) = i \h q [Scjy 1 (t) , 8-k 1 (t) ; t] , Stf (t) 



Sir 1 (t ) = i\h [8(f) 1 (t ) , 8-k 1 (t ) ; t] , Sir 1 (t ) 



(14) 



and the initial conditions 8(f) 1 (to) = 8(f>(to), Sir 1 (to) = 8n(to)- Because Ho is quadratic, the interaction picture 
operators are free fields, satisfying linear wave equations. 

It follows from Eq. (fT4)) that in evaluating Ho [8(f) 1 , Sir 1 ; t] we can take the time argument of 8(f) 1 and Sir 1 
to have any value, and in particular we can take it as to, so that HolS^ 1 ^), Sir 1 (t); t] — Ho[8<p(to), Sir(to); t] , 
but the intrinsic time-dependence of Ho still remains. The solution of Eq. (|14j) can again be written as a 
unitary transformation: 



S^ a (t) = U^fatotfMtoPo&to) , SirKt) = ^'(MoMtt^o^oCMo) 
with Uq defined by the differential equation 



(15) 



— Uo(tM) = -i H [S(f)(to),8Tr(to);t}Uo(t,to) 
at 

and the initial condition Uo(to,to) = 1 • Then from Eqs. (fT2"|) and (|16l) we have 



d_ 

dt 



U -\t,t )U(t,t ) ^-iU 1 (t,t )Hi[8(f)(to),STr(to);t]U(t,to) . 



This gives 



(16) 



U(t,to) = Uo(t,t )F(t,t ) , -F(t,to) = -iHi(t)F(t,t ) , F(t ,t Q ) = 1 

dt 

where Hi(t) is the interaction Hamiltonian in the interaction picture: 

Hi(t) = U Q (t, t Q )Hi[S(f)(to), 8n(to);t]U^(t, to) = Hi[8^(t), Sir 1 \t);t] 



(17) 



(18) 
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The solution of equations like (|17p is well known (see for example [72| ) 

F(t,t ) = Texp {^i H!(t)dt^ , 



(19) 



where T indicates that the products of His in the power series expansion of the exponential are to be time- 
ordered; that is, they are to be written from left to right in the decreasing order of time arguments. The 
solution for the fluctuations in terms of the free fields of the interaction picture is given by Eqs. (fTTj) and 
([T9"f . Then expectation values of some product A(t) of field operators are obtained as 



(A(t)) 



Texp 



Hi{i)dt 



Texp 



ffj(i) dt 



to 



(20) 



where A(t) is any 5<fi(x,t) or dir(x, t) or any product of the 5<fis and/or Sns, all at the same time t but in 
general with different space coordinates, and A 1 (t) is the same product of <5</> 7 (x,i) and/or <57r 7 (x, t). Also, 
T denotes anti-time-ordering: products of H/s in the power series expansion of the exponential are to be 
written from left to right in the increasing order of time arguments. It is more convenient to use a formula 
equivalent to Eq. (|2"0"|) : 



(At)) = E 



N 



N=0 



dt 



N 



dt 



JV-l ' ' ' 



dh 



Hi{tx\ \H I (t 2 ) 1 ---\H I {t N ) 1 A I (t) 



(21) 



In the following section, we apply this formula to calculate the bispectrum of quantum fields generated 
during inflation. 



III. NON-LINEAR COSMOLOGICAL PERTURBATIONS 



In this section, we calculate the action for higher order cosmological perturbations that is necessary to 
compute the interaction Hamiltonian at the third order in the in-in formalism. This calculation is pioneered 
by Ref. 0] and extended to general inflation models by Ref. [l8|, [III [24j ■ Here we review the derivation of 
the higher order action by following Ref. [24|, HH . 



A. Inflation models 



To set up our notation, let us first review the formalism in [73j where a general Lagrangian for the inflaton 
field is considered. The Lagrangian is of the general form 

S = 1 J d A x^-g [M 2 pl R + 2P(X, 0)] , (22) 

where (f> is the inflaton field and X = —(l/2)g fJ, "d fl cf)dv4>. The reduced Planck mass is M p i = (8ttG)-^ 2 and 
the signature of the metric is (—1,1,1,1). The energy of the inflaton field is 

E = 2XP^ X - P , (23) 

where Pjc denote the derivative with respect to X. Suppose the universe is homogeneous with a Friedmann- 
Robertson- Walker metric 



ds z = -dt z +a\t)dxi, . (24) 

Here a(t) is the scale factor and H = a/a is the Hubble parameter of the universe. The equations of motion 
of the gravitational dynamics are the Friedmann equation and the continuity equation 

ZMlH 2 = E , E = -3H(E + P) . (25) 
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It is useful to define the "speed of sound" c s as (see [z3, [75[ for a definition of the sound speed) 



and some "slow variation parameters" as in standard slow roll inflation 

H 2 M^H 2 ' 7 ~ eH ' c s H ■ ( ' 

These parameters are more general than the usual slow roll parameters (which are defined through properties 
of a flat potential, assuming canonical kinetic terms), and in general depend on derivative terms as well as 
the potential. For example, in Dirac-Born-Infeld (DBI) inflation the potential can be steep, and kinetically 
driven inflation can occur even in absence of a potential. We also note that the smallness of the parameters 
e, i], s does not imply that the rolling of inflaton is slow. When we refer to the slow-roll expansion, we 
assume that all the three slow variation parameters are small. 

The primordial power spectrum is derived for this general Lagrangian in [73| 

. . 1 P 2 1 P 2 . . 

Vd ) = 36n 2 M*c s (P + E) = 8^M^^ ' (28) 

where the expression is evaluated at the time of horizon exit at c s k = aH and Pf(fc) = k 3 P^(k) / 2ir 2 . The 
spectral index is 

n s ~l= - C y =-2e-T}-s . (29 
dm fc 

In order to have an almost scale invariant power spectrum, we need to require the 3 parameters e, ry, s to be 
very small, which we will denote simply as 0(e). We note that in inflationary models with standard kinetic 
terms the speed of sound is c s — 1. In the case of DBI inflation, the speed of sound can be very small. In the 
case of arbitrary c s , Eqs. ([28} and (f29|) for the power spectrum and its index at leading order is still valid as 
long as the variation of the sound speed is slow, namely s <§; 1. In the following we set M p i = 1. 



B. Effective action for higher order perturbations 

Now in the general setup described by the action (|22|) , we need to expand the action up to the cubic order 
in perturbations to obtain the third order interacting Hamiltonian. For this purpose,it is convenient to use 
the ADM metric formalism [76(. The ADM line element reads 

ds 2 = -N 2 dt 2 + hij [dx i + N l dt) {dx J + N ] dt) , (30) 

where N is the lapse function, N l is the shift vector and is the 3D metric. 
The action (|22|) becomes 

S = i J dtd 3 xVhN ( (3) i? + 2P) + 1 J dtd^VhN- 1 {E ij E ij - E 2 ) . (31) 

The tensor E^ is defined as 

E i} = \ (hii ~ ViNj - VjNi) , (32) 

and it is related to the extrinsic curvature by = N~ l Eij. Vj is the covariant derivative with respect to 
hij and all contra-variant indices in this section are raised with h lJ unless stated otherwise. 
The Hamiltonian and momentum constraints are respectively 

(3) i? + 2P - 2n 2 N- 2 P. x - N- 2 {E^E 1 ' 3 - E 2 ) = 0, 

Vj (iV- 1 ^') - V, (N-'E) = nN^Vi^Px, (33) 
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where 7r is defined as 

7T = - 2V 3 'Vj0. (34) 

We decompose the shift vector N 1 into scalar and intrinsic vector parts as 

Ni = N,+ d^, (35) 

where diN 1 — 0, here indices are raised with £y. 

Before we consider perturbations around our background let us count the number of degrees of freedom 
(dof) that we have. There are five scalar functions, the field 0, N, if), deth and ~ didjH, where H is a 
scalar function and deth denotes the determinant of the 3D metric. Also, there are two vector modes N l and 
hij ~ diXji where x J is an arbitrary vector. Both and x 3 satisfy a divergenceless condition and so carry 
four dof. Furthermore, we also have a transverse and traceless tensor mode 7y that contains two additional 
dof. Because our theory is invariant under change of coordinates we can eliminate some of these dof. For 
instance, a spatial reparametrization like x l — x 1 + d l e{x, t) + e % , t ->{x, i), where e and are arbitrary and 

di e \t) = ^' can ^ e cnosen so that it removes one scalar dof and one vector mode. A time reparametrization 
would eliminate another scalar dof. Constraints in the action will eliminate further two scalar dof and a 
vector mode. In the end we are left with one scalar, zero vector and one tensor modes that correspond to 
three physical propagating dof. In this review, we are primarily interested in a scalar degree of freedom. 

In order to identify this scalar degree of freedom, we need to fix a gauge. There are two commonly used 
gauges. In the next subsection we derive the higher order action in these gauges. 

C. Non-linear perturbations in the comoving gauge 

In the comoving gauge, the scalar degree of freedom is the so-called curvature perturbation £ and the 
inflaton fluctuations vanish. The 3D metric is perturbed as 

5<j) = 0, 

= a 2 e 2C h i: j , hij = + 7y + ^7ifc7i H ( 36 ) 

where deth = 1, 7^ is a tensor perturbation that we assume to be a second order quantity, i. e. 7^ = 0(C, 2 ). 
It obeys the traceless and transverse conditions 7] = d*7y = (indices are raised with S^). £ is the gauge 
invariant scalar perturbation. In (f36)) . we have ignored the first order tensor perturbations '^ygy. This 
is because any correlation function involving this tensor mode will be smaller than a correlation function 
involving only scalars, see results of [3j. 

We expand N and N l in power of the perturbation £ 

N = 1 + at + a 2 + ■ ■ ■ , (37) 

Ni = Ni (1) + A> ( (2) + ■ • • , (38) 
i> = 1>i + V>2 + • • • , (39) 

~ (™) 

where a„, Ni and ipn are of order In order to compute the effective action to order C(C 3 ), as pointed 
out in Q , in the ADM formalism one only needs to consider the perturbations of N and N l to the first order 
0(C)- This is because their perturbations at order C(C 3 ) such as a% will multiply the constraint equation 
at the zeroth order C(C°) which vanishes, and the second order perturbations such as a.2 will multiply a 
factor which vanishes by the first order solution. So the first order solution for TV and N l is enough for our 
purpose. Therefore our task is simplified. In order to expand the action (|22p to quadratic and cubic order 
in the primordial scalar perturbation £, we only need to plug in the solution for the first order perturbation 
in N and N l and do the expansion. 

Now, the strategy is to solve the constraint equations for the lapse function and shift vector in terms of 
C and then plug in the solutions in the expanded action up to third order. At first order in £, a particular 
solution for equations (|33|) is d, : 

«i = i, ^ (1) -0, Vi = -4+X. 3 2 X = a 2 4c ( 4 °) 
H H cj 
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The second order action is 



S 2 = / dtdfx 



o 3-Lc a -ae(0C) a 



The third order action is 0, 0, [HI 



S3 = / dtd 3 x 



-eaC(aC) 2 -a 3 (S + 2A)-| J + ^CC 2 



Here we defined two parameters following [IE 



H 2 e 

n 2 



X — X P t xx + P,xxx 



(41) 



(42) 

(43) 
(44) 



D. Non-linear perturbations in the uniform curvature gauge 

In this gauge, the inflaton perturbation does not vanish and the 3D metric takes the form 

<f>{x,t) = 0o + 5(j){x,t), 



1 



hij = a z h l3 , hij = Sij + jij + -Jiklj + 



(45) 



where deth = 1 and 7^ is a tensor perturbation that we assume to be a second order quantity, i.e., 7^ 
(D(5(f> 2 ). It obeys the traceless and transverse conditions 7? = 9 l 7y = (indices are raised with Sij). 
We expand N and N l in powers of the perturbation S(p(x, t) 



N = 1 + a x +a 2 H 

iV* = tv 4 (1) +iv, (2) +■•• 

^ = ipi + ip 2 H 



(46) 

(47) 
(48) 



where a„, iVj and ^>„ are of order i50 n and (j>g(t) is the background value of the field. At first order in 
a particular solution for equations (|33p is 0, : 



1 • 



A~r (!) 



2P 



The second-order action is given by 



c 2 dt 



(49) 



dtd 3 xa 3 



P.xxX { 



2X oai 



-P 



+ \ [ ^ - 0oai^ + X a 2 - a" 2 Q(cW0) 2 + ^50^ 



-3H 2 al + P$ai64> - 2 a - 2 Ha 1 d 2 tp 1 , 



(50) 
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where Xq = 4>q/2. The third-order action is obtained as 



S 3 = I dtd 3 xa 3 



Pxx I ^o^ 3 ■+ -V, 



-2 ' 1 • 

-a . 

2 



ooji (—46<f> + 5<t>aa.x8<t> — AJCqo^ 

^ 5(j){d5(j)) 2 + X ai{dS(j)) 2 - 2X Q (s'<t> - 0oai) d % 80^ 



+P.X4, (^HH 2 - foa^U + X a\54> - a~ 2 Q<ty(0ty) a + ^8^80^ 

(1 . . 3 / . 2 • - 4 

-00^ +X o ai(-2(50 + 2^ o ai50 - -^o«i 

2 „ i „ , o „ A „ / 1 • • , 



+-Pxx<^o (<5<M0 - 2(j> ax5(f>54> + 2X al8(f>J + P t x$<j> ( - Xo«i 



, 2 



+P,x ai 



2 f 



+ 0oO!i<ty - X a 2 j - aT 2 ( ^ ai (d5(f>) 2 + (s'(j> - <p Q a^j diS0ip! 



+\P l 4,4,oi 1 5(j> 2 + \p^5cj? + 3H 2 af + 2 a - 2 Ha 2 d 2 ^ 1 + ~a"V ((9 2 Vi) 2 - c^iM^i) 



(51) 



E. Relation between gauges 

The gauges used in the previous two sections are of course related by a gauge transformation. Introducing 
a new variable defined by £„ = —H8<f>/<j)Q, ( in the comoving gauge is related to 8(f) in the flat gauge as [3j 

C = Cn + /(Cn), (52) 

where 

/(C) = ^e + 4r T ti + T^\-(dO{dQ + d- 2 (d l d ] {d. l td J Q)\ +-±-\(d()(d x )-d- 2 (d l d J (d l <;d 3 x)) 



4c? c?H 4a 2 # 2 



2a 2 ff 



(53) 

On large scales where £ n becomes constant and we get 

C = Cn + ^C*. (54) 

We can show that this is nothing but the expression for £ obtained in the delta-N formalism using the 
relations 

N* = ~ = (55) 

<Po <Po <Po 

and the definition of T) in Eq. (|2?p . 

There are two ways to calculate the bispectrum of £ on large scales. One is to calculate the bispectrum 
of ( n = —H8(f)/4>o in the flat gauge and apply the delta-N formalism. It is also possible to calculate the 
bispectrum of ( in the comoving gauge. We will use both approaches in the next section. 

IV. BISPECTRUM OF CURVATURE PERTURBATION 



In this section, we first calculate the three point function for the field perturbations in the in-in formalism 
using the cubic order action obtained in the previous section. We only consider leading order terms in 
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slow-roll expansions. In the following, we consider two inflation models, k-inflation and standard slow-roll 
inflation. Then using the delta-N formalism, we derive the bispectrum of the curvature perturbation. We 
follow the calculations in Ref. 0, HH . We also discuss a method to calculate it directly in the comoving 
gauge presented in Ref. [24| . 



A. Bispectrum of quantum fields in k-inflation 



First let us consider models with non-standard kinetic terms. This is known as k-inflation. In this case, 
the leading order terms in the slow-roll expansion in the action in the flat gauge (|50[) and (1511) are given by 



dtd 3 x 



a 3 Px 



1 



1 



{dS(f>) 2 



(56) 



S3 = / dtd A x 



Pxx^aWidScj)) 2 



(57) 



The perturbations in the interacting picture are promoted to quantum operators like 

a0( T X ) = _L_ f d 3 kc50(T,k)e ,;k x , <ty(T,k) =u(r,k)a(k)+ii*(r,-k)a t (-k). (58) 
(Air)* J 

a(k) and a^(— k) are the annihilation and creation operator respectively, that satisfy the usual commutation 
relations 

[a(kO,at(k a )] =(2 7 r) 3 ^ 3 )(k 1 -k 2 ), [a(ki),o(k a )] = [ot(ki),ot(k a )] =0. (59) 
At leading order the solution for the mode functions is given by 

H 1 



u(t, k) 



(1 + ikc s r) e 



-ikc s r 



(60) 



Using the in-in formalism Ec 



the vacuum expectation value of the three point operator in the 



interaction picture is written as [3j, |7lj 

(^(i,k 1 )^(i ) k 2 )^(t,k 3 )) =-i I ^<[^(i,ki)^(t 1 ka)^(t,k J ,),fl/(t)]), 



(61) 



to 



where to is some early time during inflation when the field's vacuum fluctuation are deep inside the horizons, 
t is some time after horizon exit. If one uses conformal time, it's a good approximation to perform the 
integration from — oo to because r m —(aH) . Hi denotes the interaction Hamiltonian and it is given 
by Hi = — L3, where L3 is the Lagrangian obtained from the action (1571) . using the solution for the mode 
function and commutation relations for the creation and annihilation operators, we get 



(<ty(ki)<ty(k a )(ty(k3)) = -(27r) 3 «^(k 1 +k 2 +k 3 ) 



H 4 



where 



A 



k—inf 



3 .A k-y k*2 



E A' 3 



fcfk2_k3 

A 



2e C 2(P x ) 3 /2IT 3 =1 fc 3 



k-2 + ka k-ika 
1+ A + 2 



— mf (k 1 ,k 2 ,k 3 ), 



A 



K 2 



2 cyclic terms 



2 A \ 3/^/^2 ^3 1 — I 1 ^--v 2 1 2 / 2 7 3 ^ V~~^ 7 3 

7 \ i>J ' * 



(62) 

(63) 
(64) 
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B. Bispectrum of quantum fields in slow-roll inflation 



In slow-roll inflation with a standard kinetic term P(X) = X — V(cj>), leading-order terms in the third-order 
action (|5Tj) are given by 

dtdx 3 a 3 ( - ^5<t>dipd6<t> - ^^(U) 2 - ^^^(ddcP) 2 ^ , 

where ip was defined in Eq. (I49p and is given to leading order by 

5V = (65) 
The three point function can be calculated in the same way. We obtain 

W(k 1 )^(k 2 )^(k 3 )) = -(27r) 3 <5( 3 >(k 1+ k 2 +k 3 ) J t n3 1 g -A^Cfci.fea.fcs), (66) 

where 

Af™*^, fca, to) = e ( -± £ *? + ± £ fc * fc ' + Jc E ^ fc * I ■ ( 6? ) 



C. Bispectrum of curvature perturbation 



Now we can apply the delta-N formalism to calculate the bispectrum of the curvature perturbation. We 
define the bispectrum of the curvature perturbation as 

(C(k 1 )C(k 2 )C(k 3 )) = (2^) 7 ,5 3 (k 1 + k 2 + k 3 )(7> c ) 2 ^3 A c (h,k 2 , k 3 ) , (68) 

where Vq is given by Eq. (1281) . In k-inflation models, in the small sound speed limit and at the leading 
order in slow-roll expansion, the relation between the curvature perturbation and the field perturbation is 
simply given by £ = ~H5<fi/(j> = P^xS(j)/2e. Then the three point function for £ is given by Eq. (155)) where 
Aq = A k f mf given by Eq. (j64]). 

In standard slow- roll inflation case, the relation between £ and £„ can be written as £ = ( n + rjC^/A. Then 
the the bispectrum of the curvature perturbation is given by Eq. (|68[) where ^4^ is given by 




D. Computation in comoving gauge 



It is also possible to calculate the bispectrum of the curvature perturbation in comoving gauge and this 
gives a very useful consistency check. Here we follow Ref. [24j and see how this works. 

In fact the cubic effective action in (|57l) looks like order O(e ) in the slow variation parameters while in 
the previous section, we find that the bispectrum is suppressed by slow-roll parameters in slow-roll inflation. 
In slow-roll inflation, as emphasized and demonstrated in Ref. [3j, one can perform a lot of integrations by 
parts and cancel terms of order <D(e°) and 0(e). The resulting cubic action is actually of leading order 0(e 2 ) 
in slow roll parameters. A similar analysis can be performed for the general Lagrangian in Ref. [l8j |. Except 
for terms that are proportional to 1 — c 2 s or A, the rest of the terms can be cancelled to the second order 
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0(e 2 ) and the cubic order action Eq. (l4"2l can be rewritten as 



5, = 



dtd 3 
c 



-(^-i + 2A)^ + f(e. 



a e a 
2&~dt \ . 



( 2 ( + ^(d()(dx)d 2 X - 



4a 



(d 2 C)(dx) 2 + 2f(0 



SL 



St 



where x is defined in Eq. (|40)) and in the last term 



SL 



(70) 



(71) 



Here d~ 2 is the inverse Laplacian, 8L/SC\\ is the variation of the quadratic action with respect to the 
perturbation £, therefore the last term which is proportional to SL/5^\i can be absorbed by a field redefinition 
of £. It can be easily shown that the field redefinition that absorbs this term is 



C -► Cn + /(<„), 



(72) 



where /(C) is given by Eq. (|53[) . This is nothing but the relation between £ and £„ obtained from the guage 
transformation between the flat gauge and comoving guage. One then computes the vacuum expectation 
value of the three point function in the interaction picture in the same way. We get Eq. (|68l) with 



A 



- 1 - 



2^^\ 3/^2 ^2 ^3 



2K 3 



1 



- 1 



\ i>3 




(73) 



In k-inflation, the first two terms arc the leading order contributions in slow-roll expansions. The remaining 
terms are 0(e). Note that one should take into account 0(e) corrections from the leading order contribution 
in order to obtain a full expression up to 0(e) in k-inflation. In standard slow- roll inflation, the first two 
terms and the last term vanishes and c„ = 1. 



V. SHAPES OF BISPECTRUM 



In this section, we compare the prediction of the bispectrum in k-inflation with the local-type non- 
Gaussianity. The discussions in this section are based on Refs. (78l. [79j. 

The bispectrum in the local-type non-Gaussianity is often characterized by 

C = & + (74) 

where Q n obeys Gaussian statistics. Originally the /j^£ ai parameter was introduce to parametrise a non- 
linearity in the curvature perturbation $ in the Longitudinal gauge which is related to C as $ = (3/5)£. 
Note that Ref. [3| uses a different sign convention for /j^£ a/ from WMAP papers (see 80] for the latest 
result). Here we follow the definition used in WMAP papers. Eq. ([71)1 is nothing more than the expression 
for £ in the delta-N formalism. In slow-roll inflation /j^£ a ' is 0(e) but models like curvaton predicts /^£ o/ 
larger than one. The bispectrum of curvature perturbation is given by 



(C(k 1 )C(k 2 )C(k 3 )) = (2^) 3 <5 3 (k! + k 2 + k 3 )(V c ) 2 F(k u k 2 , fc 3 ), 



(75) 
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where 

F local {k 1 MM) = {^) i (yf^T l ) (^3+^ + ^3)- (76) 

Eq. (|74|) describes (at leading order) the most generic form of non-Gaussianity which is local in real space. 
This form is therefore expected for models where non-linearities develop outside the horizon. This happens 
for all the models in which the fluctuations of an additional light field, different from the inflaton, contribute 
to the curvature perturbations we observe. In this case non-linearities come from the evolution of this field 
outside the horizon and from the conversion mechanism which transforms the fluctuations of this field into 
the curvature perturbations. Both these sources of non-linearity give a non-Gaussianity of the form (I74[) 
because they occur outside the horizon. Examples of this general scenario are the curvaton models [81|, 
models with fluctuations in the reheating efficiency [82[ and multi- field inflationary models [f| . 

Being local in position space, Eq. (|74|) describes correlation among Fourier modes of very different k. It 
is instructive to take the limit in which one of the modes becomes of very long wavelength [3], k 3 0, 
which implies, due to momentum conservation, that the other two fc's become equal and opposite. The long 
wavelength mode £jg freezes out much before the others and behaves as a background for their evolution. 
In this limit Fi oca i is proportional to the power spectrum of the short and long wavelength modes 



Flocal OC -3-p . (77) 

«, 3 h, x 



This means that the short wavelength 2-point function (fr ) depends linearly on the long wavelength 
mode C^ 3 

(Cs.CfcCs,) <* (Cs 3 C_g 3 >^<Cs 1 C_g 1 > • (78) 

From this point of view we expect that any bispectra will reduce to the local shape ([75|) in the degenerate 
limit we considered if the derivative with respect to the long wavelength mode does not vanish. 

In standard single field slow- roll inflation, as pointed out in [3j, different points along the background 
wave are equivalent to shift in time along the inflaton trajectory, so that the derivative with respect to the 
background wave is proportional to the tilt of the scalar spectrum. This can be explicitly checked in the full 
expression of the 3-point function (Eqs. (p?) and 



(27r) 4 

F s tand(kl, k 2 , h) = -Q-p- 

In the limit k 3 — > Eq. (|79)l goes as 




(79) 



F stand (k 3 -> 0) oc (2e + 77) -3-3 = -(n. - 1 W . (80) 

3 1 3 1 

As expected the tilt in the spectrum n s fixes the degenerate limit of the 3-point function. Note however 
that expression (|79p is not of the local form (|76p but contains contributions which are important for non- 
degenerate triangles. If we compare expression (1761) and (|79[) and neglect the different shape dependence, we 
see that standard single-field inflation predicts /j^£ ai of order of the slow-roll parameters. 

We have seen that the degenerate limit k 3 —> describes the effect of a slowly- varying long-wavelength 
perturbation on the 2-point function of short wavelength modes. In many models, the correlation is much 
weaker in this limit than in the local model (|76j) . Physically this means that the correlation is among modes 
with comparable wavelength which go out of the horizon nearly at the same time. In this case the 3-point 
function in the degenerate limit is suppressed by powers of £3 with respect to the behaviour of Eq. (|77[), 
We have correlation among modes of comparable wavelength in all models in which the non-Gaussianity 
is generated by derivative interactions: these interactions become exponentially irrelevant when the modes 
go out of the horizon because both time and spatial derivatives become small, so that all the correlation is 
among modes freezing almost at the same time. 
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k, k, 



FIG. 1: Left: Plot of the function F(l, k 2 /ki, k 3 /k 1 )(k 2 /k- l ) 2 (k 3 /k 1 ) 2 predicted by the DBI models QjJ. Right: 
Difference between the above plot and the analogous one (top of fig. [5J) for the factorizable equilateral shape used in 
the analysis. From [7S|. 



K-inflation is a typical example for these type of models. The three point function is obtained in Eq. (|73[) 
and given by 



F(k u k 2 ,h) - !2;n 



1 2A\ 3fcffc|fcJ ( 1 



2 £ J 2K 3 ' V ^ 




In a model of inflation based on the DBI action, 



P(cf>,X) = -/(</>)- Vl- ZXftt) + f{^r 1 + Vty), (82) 

A is given by 

A=|fi-lV (83) 



Thus the first term vanishes in Eq. (|5Tj) . Unfortunately, a function F is not factorizable, so it is not easy to 
perform an optimal analysis using CMB observations. However, it is a very good approximation to take a 
factorizable shape function F which is close to Eq. ([8lj) and perform the analysis for this shape. In the limit 
ki — !• with k 2 and k% fixed, all the equilateral functions diverge as k^ 1 [78| (while the local form eq. ((75)) 
goes as k^ 3 ). The factorizable function that satisfies this condition is given by 

A 4 ( 9 feg uil\ /111 2 1 



F(ki, k 2 , fe 3 ) = (2tt) -/A 1 -Tars - I3T3 - I3T3 - 757512 + T-riTs + (5 Perm.) 1 , (84) 

\^_LU y \^ rii^n^ 13 2 3 12 3 1 2 3 / 

where the permutations act only on the last term in parentheses. In figurc[T]we study the equilateral function 
predicted in DBI inflation (|8Tj) . In the second part of the figure we show the difference between this function 
and the factorizable one used in our analysis. We see that the relative difference is quite small. The same 
remains true for other equilateral shapes (see [78| for the analogous plots for other models). 

In figure we compare this function with the local shape. The dependence of both functions under a 
common rescaling of all fc's is fixed to be oc fc~ 6 by scale invariance, so that we can factor out k^ 6 for 
example. Everything will now depend only on the ratios k 2 /ki and k^/ki, which fix the shape of the triangle 
in momentum space. For each shape we plot F(l, k 2 /ki, k$ / 'ki)(k 2 / 'ki) 2 '(k^ / 'ki) 2 this is the relevant quantity 
if we are interested in the relative importance of different triangular shapes. The square of this function 
gives the signal to noise contribution of a particular shape in momentum space [78| . We see that for the 
function (|54"|) , the signal to noise is concentrated on equilateral configurations, while squeezed triangles with 
one side much smaller than the others are the most relevant for the local shape. 
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Factorizable eq. shape 1 *2 Loca 1 shape 1 




FIG. 2: Plot of the function F(l, k%/k\, kz/ki){k2/ki) 2 (k3/k{) 2 for the equilateral shape used in the analysis (left) 
and for the local shape (right). The functions are both normalized to unity for equilateral configurations = = 1. 

Since F(ki, k2, k^) is symmetric in its three arguments, it is sufficient to specify it for ki > k^ > £13, so ^ < < 1 
above. Moreover, the triangle inequality says that no side can be longer than the sum of the other two, so we only 
plot F in the triangular region 1 — j^-<j^<j^<l above, setting it to zero elsewhere. From [7^ ]. 



VI. CONCLUSION 



In this review, the bispectrum of curvature perturbation is calculated using the in-in formalsim and the 
delta-N formalism. There are two distinct contributions to the bispectrum. One is coming from a non-linear 
relation between the curvature perturbation £ and quantum fluctuations of a field at the horizon crossing. 
In this case the non-linearities come from the evolution of this field outside the horizon. Being local in a 
position space, the shape of the bispectrum is highly non-local in a Fourie space having a maximum signal for 
the squeezed configuration k% -C k\ , hi . The other contribution is coming from the bispectrum of quantum 
fields generated under horizon scales. In models in which the non-Gaussianity is generated by derivative 
interactions such as DBI inflation and k-inflation models, we have correlation among modes of comparable 
wavelength in all models and these interactions become exponentially irrelevant when the modes go out of 
the horizon because both time and spatial derivatives become small, so that all the correlation is among 
modes freezing almost at the same time. Then the bispectrum has a peak at the equilateral configuration 

h ~ &2 ~ fa- 

In order to put constraints on the bispectrum from CMB observations, it is necessary to construct an 
estimator that uses a model prediction for the bispectrum as an template. For the local-type non-Gaussianity 
and the equilateral non-Gaussianity, the constraints obtained in WMAP 7-year results are [80j 

- 10 < f l £i al < 74, -214 < f*f x < 266. (85) 

at 95% confidence level. Recently, it has been found that Large Scale Structure (LSS) can give a similar 
level of constraints on the local type non-Gaussianity from the scale-dependent bias effects on the halo 
power spectrum [H| while this effect is absent in the equilateral non-Gaussianity [13, HE]- The constraint 
on f§f l is obtained from SDSS as -29 < < 70 [sj| and combining it to WMAP 7-year results, we 

get — 5 < < 59 [80]. Thus currents observations are consistent with Gaussian primordial curvature 

perturbations but future experiments such as Planck will give much tighter constraints ~ 0(1) and 

we may be able to detect a deviation from Gaussianity which has a huge impact on early universe models. 
There are a lot issues that are not covered by this review. We will mention some of the issues here: 

• Trispectrum 

In this review, we concentrated on the leading order non-Gaussianity, i.e. the bispectrum but it 
has been recognized that the trispectrum could give a useful information to distinguish between many 
possible models that predict large non-Gaussianity. For the local non-Gaussianity, we can easily extend 
the model by expanding £ up to the third order £ = C„ + + 9NL(n/&- The trispetrum is 

characterized by two parameters tnl = f^L an( i 9nl [17l.l28l.l77l]. The constraints on these parameters 
are rather weak, -631 < < 717 and -3.80 < sjvl/IO 6 < 3.88 from WMAP 5-year results 



15 



at 95% confidence level [83] ■ But again the Planck will improve these significantly [88| . The full 
trispectrum in DBI inflation at the leading order in small sound speed limit has been obtained [25], [89l — 
[9lj . Unlike the bispectrum, there are still two degrees of freedom even for the equilateral configurations 
k% ~ &2 ~ ^3 ~ ^4 and also the form of the trispectrum is too complicated to be used for the estimator. 
It is necessary to develop approximations for the shape of the trispectrum in DBI inflation as is done 
for the bispectrum. 

Multi-field inflation 

In single field inflation models with a standard kinetic term, the resulting non-Gaussianity is small 
suppressed by slow-roll parameters. However, in multi-field models, it is possible to have large local 
non-Gaussianity due to non-linear dynamics of fields outside horizon. The dclta-N formalism is easy to 
be extended to multi-field models and there have been extensive study of non-Gaussianity in multi-field 
models using the delta-N formalism (see for example [92| and references therein). Multi field effects are 
also important in DBI inflation. In DBI inflation, fluctuations along the entropy directions of the fields 
that are orthogonal to the field trajectory have the same sound speed as the adiabatic fluctuations 
along the field trajectory 30-321 [93j . If the trajectory makes a turn in a field space, this converts 
the entropy perturbations to the curvature perturbations. Although the bispectrum is enhanced by 
this conversion, the enhancement of the power spectrum is stronger and f^i? becomes smaller in 



multi-field models, which help ease stringent constraints on DBI inflation models in string theory [31 1 



It has been shown that the trispectrum is enhanced for a given 

f e NL l [Hi 1951 ■ K is also possible 
that the multi-field effects modify the bip sectrum for quantum field at the horizon crossing. In the 
so-called quasi-single inflation models [961 ] . the entropy perturbations develops large non-Gaussianity. 
The conversion of entropy perturbations to the curvature perturbation can happen near the horizon 
crossing and during this transition the shape of the bispectrum can be modified in a non-trivial way. 

There are many other possibilities to get large non-Gaussianity of quantum fields such as a feature in 
inflaton potentials (2(| [97]] • All these models predict distinct shapes of the bispectrum and trispectrum. 
In the future, we may be able to exploit CMB and LSS data to distinguish between many possible early 
universe models via non-Gaussianity. 
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